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Abstract The magnetic fields in the solar atmosphere 
structure the plasma, store free magnetic energy and 
produce a wide variety of active solar phenomena, like 
flare and coronal mass ejections(CMEs). The distri- 
bution and strength of magnetic fields are routinely 
measured in the solar surface( photosphere). There- 
fore, there is considerable interest in accurately model- 
ing the 3D structure of the coronal magnetic field using 
photospheric vector magnetograms. Knowledge of the 
3D structure of magnetic field lines also help us to in- 
terpret other coronal observations, e.g., EUV images 
of the radiating coronal plasma. Nonlinear force-free 
field (NLFFF) models are thought to be viable tools 
for those task. Usually those models use Cartesian 
geometry. However, the spherical nature of the solar 
surface cannot be neglected when the field of view is 
large. In this work, we model the coronal magnetic 
field above multiple active regions using NLFFF extrap- 
olation code using vector magnetograph data from the 
Synoptic Optical Long-term Investigations of the Sun 
survey (SOLIS)/ Vector Spectromagnetograph (VSM) 
as a boundary conditions. We compare projections of 
the resulting magnetic field lines solutions with their 
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respective coronal EUV-images from the Atmospheric 
Imaging Assembly (SDO/AIA) observed on October 11, 
2011 and November 13, 2012. This study has found 
that the NLFFF model in spherical geometry recon- 
structs the magnetic configurations for several active 
regions which agrees with observations. During Octo- 
ber 11, 2011 observation, there are substantial number 
of trans-equatorial loops carrying electric current. 

Keywords Active Regions, Magnetic Fields; Active 
Regions, Models; Magnetic fields, Corona; Magnetic 
fields. Photosphere; Magnetic fields. Models 



1 Introduction 

The geometry and dynamics of the solar corona are de- 
termined by the evolving magnetic field at the Suns 
surface (photosphere). Thus, magnetic fields are be- 
lieved to play a dominant role for active phenomena 
carried out in the solar corona. In order to study solar 
eruptive phenomena, we have to understand how mag- 
netic energy is stored in the pre-eruptive corona. Hence, 
the three dimensional (3D) structure of magnetic fields 
and electric currents in the pre-eruptive corona and 
the amount of free energy stored in the field have to 



determined (Schrijver and Title 2005 Jiang and Feng 



20121. Unfortunately, reliable magnetic field measure- 



ments are still restricted to the level of the photosphere, 
where the inverse Zeeman effect in Fraunhofer lines is 
observable. Even if the measurement of magnetic fields 
in the solar corona has considerably improved in recent 



decades (Lin et al. 2000 Liu and Lin 2008), further 



developments are needed before accurate data are rou- 
tinely available. 

As an alternative to measurements in the solar atmo- 
sphere, magnetic field extrapolation have been imple- 
mented to calculate the field from the measured photo- 
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spheric field using assumptions that the fields are force- 
free. In this model assumption, the corona magnetic 
forces are dominant that all non-magnetic forces like 
plasma pressure gradient and gravity can be neglected 
in the lowest order (Gary 2001 1. This implies that, if 



there are currents, they must be aligned with the mag- 
netic field. To describe the equilibrium structure of the 
static coronal magnetic field, the force-free assumption 
is appropriate: 



(V X B) = aB 



V • B = 



subject to the boundary condition 



B = Bn 



on photosphere 



(1) 



(2) 



(3) 



where B is the magnetic field and Bobs is measured vec- 
tor field on the photosphere. Equation ([T]) states that 
the currents are co-aligned with magnetic fields, where 
a is the torsion function that represents the propor- 
tionality between the electric current density and the 
magnetic field. Equation ^ describes the absence of 
magnetic monopoles. 

Parameter a of Equation ([T]) can be a function of po- 
sition, but the combination of Equation ([T]) and Equa- 
tion ^ (B • Va — 0) requires that a be constant along 
a given field line, a = yields the current-free poten- 
tial field. If a is constant everywhere in the volume 
under consideration, the field is called linear force-free 
field (LFFF), otherwise it is nonlinear force- free field 



(NLFFF) (Wiegelmann"2004'; 'Wheatland and Regnier 


2009 


Wheatland and Leka||2011| |Wiegelmann et al. 


2012 


Wiegelmann and Sakurai |2012 Malanushenko 


et al. 


2012|. Potential and linear force- free fields can 



be used as a first step to model the general structure of 
magnetic fields in the solar corona, the use of nonlinear 
force free fields is essential to understand eruptive phe- 
nomena. Practically there are both observational and 
theoretical reasons which suggest that the pre-eruptive 
magnetic fields are non-linear force-free fields. For de- 
tails of those models we direct the readers to |Wiegel-| 
mann and Sakurai (2012). 



Most of the NLFFF modeling codes are implemented 
in the Cartesian geometry. Therefore those codes are 
not well suited for larger domains, since the spherical 
nature of the solar surface cannot be neglected when the 



field of view is large. DeResa et al. ( 2009 1 has studied 



that different NLFFF models have different field line 
configurations in the coronal volume with widely vary- 
ing estimates of the magnetic free energy. The study 
suggested that the main reasons for that problem are 



(1) non- force free nature of the field within the pho- 
tosphere, (2) the uncertainties on vector-field measure- 
ments, particularly on the transverse component, and 
(3) the needs of large domain that can capture the con- 
nections of an active region to its surroundings to ac- 
commodate more electric currents associated with the 
full active regions of interest. Solar Dynamics Obser- 
vatory (SDO) mission has made repeated observations 
of large, almost global scale events in which large scale 
connection between active regions may play fundamen- 
tal role. Therefore, this needs motivate us to implement 
a NLFFF procedure in spherical geometry to large scale 
boundary data from the Helioseismic and Magnetic Im- 
ager (HMI) on board SDO and Synoptic Optical Long- 
term Investigations of the Sun survey (SOLIS)/ Vec- 
tor Spectromagnetograph ( VSM) ([Wiegelmann 



Tadesse et all [20091 |20TT| [2012a|b[ [Quo et al. 
Tadesse et al.|[2013[ ) 



2007 



2012 



Coronal magnetic fields manifest themselves in X- 
rays and EUV images in the shapes of coronal loops. 
Those loops provide an additional constraint which is 
not currently being used due to the mathematical com- 
plications of incorporating such input into numerical 
models (Malanushenko et al. 2012 1. The 3D structure 



of magnetic loops from NLFFF models help us for mod- 
eling of plasma loops, and understanding coronal heat- 
ing and plasma flows along the loops. In this work, 
we apply our spherical NLFFF procedure to a group 
of active regions observed on October 11, 2011 and 
November 13, 2012 using SOLIS/VSM data. During 
the former observation, there were six active regions 
(ARs 11610, 11611, 11612, 11613, 11614 and 11615) 
and there were three vertically aligned active regions in 



later one. In our previous work Tadesse et al. (2012b), 
we have studied the connectivity between three neigh- 
boring active regions which were horizontally aligned. 
Here, we compare the extrapolated magnetic loops with 
their respective extreme ultraviolet (EUV) observations 
by the Atmospheric Imaging Assembly (AIA) on board 
SDO. 



2 Method 

This study requires extrapolating the three-dimensional 
NLFFF coronal fields from the photospheric boundary 
data. The photospheric vector magnetograms, obtained 
by the Synoptic Optical Long-term Investigations of 
the Sun survey (SOLIS)/ Vector Spectromagnetograph 
(VSM) are used as the boundary conditions. Mean- 
while, those measured data are inconsistent with the 
above force-free assumption. Therefore, one has to ap- 
ply some transformations to these data before nonlin- 
ear force-free extrapolation codes can be applied. This 



3 



procedure is known as preprocessing. This preprocess- 
ing scheme removes forces and torques from the bound- 
ary and approximates the photospheric magnetic field 
to the low plasma-/3 force-free chromosphere. For de- 
tailed descriptions of the method we direct the readers 
to Wiegelmann et al. ( 2006 1 for Cartesian geometry and 



Tadesse et al. ( 2009 ) for spherical one. 

We solve the force- free equations ([T])(with the form 
V X B) X B = 0) and (121) using optimization principle 



(Wheatland et al.|2000 


Wiegelmann|2004|) in spherical 


geometry ( 


Wiegelmann|2007 


Tadesse et al.|2009 2013 



2012c). The method minimizes a joint measure 
of the normalized Lorentz forces and the divergence of 
the field throughout the volume of interest, V. Opti- 
mization procedure is spherical geometry is given by: 



Lui — Lf + Ld + vLphoto (4) 

Lf^ ( w/(r, 6*, 0)5-2 1 (V X B) X BjV^ sin 6'(ird6'd(/) 
Jv 

Ld= I ujd{r,e,(j))\\7 -Bl^r^smedrdedcj) 
Jv 



-'photo 



I Bd.W(0, 

Js 



■ Bd r sin 



where B^; = (B — Bobs). Lf and measure how well 
the force-free Eqs. (jlj and divergence-free ^ condi- 
tions are fulfilled, respectively, and both ^, (j)) and 
ijJd(r, 9, 4>) are weighting functions. For detailed descrip- 
tion of the use of weighting functions, we direct the 
readers to Wiegelmann ( 2004 1 . The term Lphoto is sur- 



face integral over the photosphere which allows us to 
relax the field on the photosphere towards force-free 
solution without to much deviation from the original 
surface field data, 'Robs- In this integral, 'W{9,(l)) is a 
space-dependent diagonal matrix which gives different 
weights for observed surface field components depend- 
ing on its relative accuracy in measurement. 

The full inversion of SOLIS data is performed in 



the framework of Milne-Eddington model (ME)(jUnno 
1956). The inversion is made only for pixels whose 
polarization signal is above a selected threshold. Pix- 
els with polarization below threshold are left undeter- 
mined. These data gaps represent a major difficulty for 
existing magnetic field extrapolation schemes. There- 
fore, those missing data is considered most inaccurate 
and are taken account of by setting W(0, (fi) to zero in 
all elements of the matrix. For the detailed description 



of the method we direct the readers to Wiegelmann and 



Inhester (2010) and Tadesse et al. (2011) for Cartesian 



3 Results 

In this study, we have used our spherical NLFFF op- 
timization procedure to a group of active regions ob- 
served by SOLIS on October 11, 2011 and November 
13, 2012. During the former observation, there were six 
active regions (ARs 11610, 11611, 11612, 11613, 11614 
and 11615) and there were three vertically aligned ac- 
tive regions(ARs 11314, 11316, and 11319) in later one 
(see Fig. [T^ and b). In order to accommodate the con- 
nectivity between those ARs and their surroundings, we 
adopt a non uniform spherical grid r, 0, and in the di- 
rection of radius, latitude, and longitude, respectively. 
Before applying our spherical NLFFF, we preprocessed 
the two photospheric vector field data sets using our 
spherical preprocessing routine. 

There are no vector magnetic field measurements 
for the side and top boundaries of a localized domain. 
Therefore, we have to make assumptions about these 
fields before performing a NLFFF extrapolation. We 
assumed the lateral and upper boundaries of the com- 
putational domain as current-free. In order to initial- 
ize our NLFFF code, we calculated potential field from 
those data sets using their respective preprocessed ra- 
dial field components (B^) using spherical harmonic ex- 
pansion . The computational box is a wedge-shaped 
volume V with six boundary surfaces (four lateral side 
boundaries, top and photospheric boundaries). This 
box enables us to study the connectivities between ac- 
tive regions and their surroundings for the large field- 
of- views. 

For computing the NLFFF field solutions for the 
two data sets, we use the preprocessed surface vector 
boundary fields with radial, longitudinal, and latitu- 
dinal components. We minimize the functional L^^ of 
Eq. (|4|). We implement the new term Lphoto in Eg- ^ 



to treat those data gaps in SOLIS (Wiegelmann and 



Inhester 2010 Tadesse et al. 2011). For those pixels 



and spherical geometry, respectively. 



for which "Robs was successfully inverted, we allow de- 
viations between the model field B and the input fields 
observed Bq^s surface field using Eq. (|4|, so that the 
model field can be iterated closer to a force-free solu- 
tion even if the observations are inconsistent. We have 
used Langrangian multiplier oiv — 0.001 to control the 
speed with which the lower boundary is injected during 
the NLFFF extrapolation. 

In order to compare our reconstructions with obser- 
vation, we plot the selected field lines of the NLFFF 
solutions for the two data sets and we overlay the field 
lines with corresponding AIA 193 A image (see Fig- 
ure ^ and d). From Figure ^ and d, one can see 
that the field lines of NLFFF model solutions are recon- 
structed in such away that they fit with the observation. 
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Fig. 1 : a) Full-disk SOLIS/VSM magnetogram of October 11, 2011 at 18:45UT, b) Full-disk AIA 171 A image 
of October 11, 2011 at 18:45UT, c) Full-disk SOLIS/VSM magnetogram of November 13, 2012 at 20:11UT, d) 
Full-disk AIA 171 A image of of November 13, 2012 at 20:10UT. The blue boxes show the active regions observed. 
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{c) NLFFF soliuioii of October (d) NLFFF eohuion of November 13, 20l2 at 2i);llUT 



Fig. 2 : a) SDO/AIA 171 A image on October 11, 2011 at 18:45UT, b) SDO/AIA 171 A image on November 13, 
2012 at 20:10UT, c) NLFFF model field lines of October 11, 2011 at 18:45UT overlaid on AIA 193A image and 
d) NLFFF model field lines of November 13, 2012 at 20:11UT overlaid on AIA 193A image. Green and red lines 
represent open and closed magnetic field lines. 
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There are spatial correspondence between the overaU 
shape of the NLFFF field lines and the EUV loops for 
both observational dates. Those qualitative compar- 
isons between NLFFF model magnetic field lines and 
the observed EUV loops of AIA images indicate that 
the NLFFF model provides a more consistent field for 
large field-of- views. Even if the qualitative nature is 
not ideal, it remains the best option in the absence of 
a more reliable quantitative comparison. 

Once we have determined the 3D structure of the coro- 
nal magnetic field using the sophisticated numerical 
modeling, we can use those field solutions to calcu- 
late different quantities in the solar corona. We esti- 
mate the free magnetic energy to be the difference be- 
tween the extrapolated NLFFF and the potential field 
with the same normal boundary conditions in the pho- 



tosphere Regnier and Priest (2007); Thalmann et al. 



(2008). We therefore estimate the upper limit to the 
free magnetic energy associated with coronal currents 
of the form 



1 

8^ 



Iff 



B 



pot 



^sin0drd9d(j), 



(5) 



Our result for the estimation of free-magnetic energy 
in Table [l] shows that the NLFFF solutions have 2.13% 
and 6.28% more energy than the corresponding poten- 
tial field solutions obtained from data observed on Oc- 
tober 11, 2011 and November 13, 2012, respectively. 

we have 



In our previous work Tadesse et al. (2012c 



studied the connectivity between two active regions one 
in the northern hemisphere and the other in the south. 
In this study we have used the same method to calcu- 
late the percentage of the total electric current shared 
between the AR 11319 in the northern hemisphere and 
AR 11316 in the southern hemisphere observed on Oc- 
tober 11, 2011 at 18:45UT (see Fig. [3|. In order to 
quantify the percentage share in the electric current, 
we first identified those field lines carrying electric cur- 
rents and emanating from AR 11319 and ending into 
AR 11316. The ratio of total electric current density 
flux due to those electric current carrying field lines con- 
necting the two ARs to the total electric current density 
flux due to all field lines with current emanating from 
AR 11319 gives us the percentage share in the electric 

Table 1 : The magnetic energy associated with ex- 
trapolated 3D NLFFF field configurations from SOLIS 
data sets of the two dates. 



current between the two active regions. For this case 
study, we found that 23.31% of positive/negative polar- 
ity of AR 11319 in the northern hemisphere is connected 
to positive/negative polarity of AR 11316 in southern 
hemisphere. There are more electric current connec- 
tivity between these ARs as they are closer to each 
other even if they are located in different hemispheres. 
Therefore, during this particular observation, there are 
substantial number of trans-equatorial loops carrying 
electric current. 




Date 


^total(10''erg) 


i?f^cJl033erg) 


October 11, 2011 


5.63 


0.12 


November 13, 2012 


9.72 


0.61 



Fig. 3 : Selected field lines of NLFFF model solutions 
of October 11, 2011 at 18:45UT rotated to limb. 



4 Conclusion and outlook 

Coronal magnetic field modeling codes in Cartesian ge- 
ometry are not well suited for larger domains, since 
the spherical nature of the solar surface cannot be ne- 
glected when the field of view is large. In this study, we 
have investigated the 3D coronal magnetic field associ- 
ated with two data sets observed on October 11, 2011 
and November 13, 2012 by analyzing SOLIS/VSM data 
using NLFFF model. During the former observation. 



there were six active regions (ARs 11610, 11611, 11612, 
11613, 11614 and 11615) and there were three vertically 
aligned active regions in later one. We have used our 
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spherical NLFFF codes to compute the magnetic field 
solutions for large field of views with six active regions 
and three vertically aligned ones for the first time. 

For the comparison with observations, we overlayed 
the NLFFF model field solutions for the two data sets 
with their respective EUV loops from AIA. The study 
indicates that there are significant agreements between 
EUV loops and NLFFF model solutions. There are 
spatial correspondence between the overall shape of the 
NLFFF field lines and the EUV loops for both obser- 
vational dates. Those qualitative comparisons between 
NLFFF model magnetic field lines and the observed 
EUV loops of AIA images indicate that the NLFFF 
model provides a more consistent field for large field- 
of-views. Even if the qualitative nature is not ideal, it 
remains the best option in the absence of a more reliable 
quantitative comparison. 

Today, Solar Dynamics Observatory (SDO) mission 
has been repeated observations of large, almost global 
scale events in which large scale connection between ac- 
tive regions may play fundamental role. Therefore, it 
is useful to implement a NLFFF procedure in spherical 
geometry for use when large scale boundary data are 
available, such as from the Helioseismic and Magnetic 
Imager (HMI) on board SDO. In order to handle large 
resolution data from HMI, one can adapt the multi-grid 
technique to the spherical code to improve its perfor- 
mance. 
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